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Abstract 

We study the relativistic quark-antiquark system embedded in 
magnetic field (MF). The Hamiltonian containing confinement, color 
Coulomb and spin-spin interaction is derived. We analytically follow 
the evolution of the lowest neutral meson state as a function of MF 
strength. Calculating the color Coulomb energy (Vcoui) we have ob- 
served the unbounded negative (at least in the limit of large Nc) contri- 
bution at large MF which makes the mass negative for eB > eB^^^ . 
We display the vr'^ and masses as functions of MF in comparison 
with recent lattice data. 

Pacs: 



1 Introduction 



During the last years we have witnessed an impressive progress of the funda- 
mental physics in ultra-intense magnetic field (MF) reaching the strength up 
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to eB ~ lO^^G ~ P]. Until recently magnetars [2] were the only physical 
objects, where such, or somewhat weaker MF could be realized. Now MF of 
the above strength and even stronger is within reach in peripheral heavy ion 
collisions at RHIC and LHC [3j . High intensity lasers is another perspective 
tool to achieve MF beyond the Schwinger limit On the theoretical side 
a striking progress has been achieved along several lines. It is beyond our 
scope to discuss these works or even present a list of corresponding refer- 
ences. We mention only two lines of research which have a certain overlap 
with our work. The first one [5l |6] is the behavior of the hydrogen atom and 
positronium in very strong MF. The second one [7j is the conjecture of the 
vacuum reconstruction due to vector meson condensation in large MF. The 
relation between the above studies and our work will be clarified in what 
follows. 

Our goal is to study from the first principles the spectrum of a meson 
composed of quark-antiquark embedded in MF. Use will be made of Fock- 
Feynman- Schwinger representation (see |S] for review and references) of the 
quark Green's function with strong (QCD) interaction and MF included. An 
alternative approach could have been Bethe-Salpeter type formalism. How- 
ever, for the confinement originating from the area law of the Wilson loop, 
the use of the gluon propagator is inadequate. Numerous attempts in this 
direction failed because of gauge dependence and the vector character of the 
gluon propagator, while confinement is scalar and gauge invariant. Therefore 
it is sensible to use the path integral technique for QCD+QED Green's func- 
tions. This method in combination with the einbein technique (the method 
of effective masses) [9] enables one to construct explicit expressions for meson 
Hamiltonians without MF [TU]. In this way spectra of light-light, light-heavy 
and heavy-heavy mesons were computed with a good accuracy, using the 
string tension a, strong coupling constant as and quark current masses as 
an input [llj.tl2]. 

In what follows we expand this technique to incorporate the effects of 
MF on mesons. The latter contains: 1) direct influence of MF on quark and 
antiquark, and 2) the influence on gluonic flelds, e.g., on as, gluon propagator 
and on the gluon fleld correlators determining the string tension a. However, 
since MF acts on charged objects, its influence on the gluonic degrees of 
freedom enters only via (A^c)"^, k = 1, 2, ... In what follows the corrections 
of the second type will be neglected. 3) As will be discussed elsewhere, MF 
also changes quark condensate {qq) and quark decay constants /,r etc., and 
in this way strongly influences chiral dynamics. 
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The decisive step in our relativistic formalism is the implementation of 
the pseudomomentum notion and cm. factorization in MF, suggested in the 
nonrelativistic case in p3J for neutral two particle systems. 

The plan of the paper is the following. Section 2 contains a brief peda- 
gogical reminder of how the two-body problem in MF is solved in quantum 
mechanics. The central point here is the integral of motion ( "pseudomomen- 
tum") which allows the separation of the center of mass. Here we also show 
how to diagonalize the spin-dependent interaction. In section 3 we formulate 
the path integral for quark-antiquark system with QCD+QED interaction. 
Then from Green's function the relativistic Hamiltonian is obtained. Sec- 
tion 4 is devoted to the treatment of confining and color Coulomb terms, 
we demonstrate the unboundedness of qq spectrum due to the latter. This 
phenomenon can be called the "magnetic QCD collapse", which occurs in 
the large Nc limit: Uf/Nc — )■ 0. Here we also present the derivation of the 
eigenvalue equations for the relativistic Coulomb problem. In section 5 we 
discuss the spectrum of the system focusing on the regime of ultra-strong 
MF. Section 6 contains the discussion of the results, comparison with lattice 
calculations, drawing further perspectives and intersections of our results 
with those of other authors O El [7j. 

2 Pseudomomentum and Wavefunction Fac- 
torization 

The total momentum of mutually interacting particles with translation 
invariant interaction is a constant of motion and the center of mass mo- 
tions can be separated in Schroedinger equation. It was shown [13j that a 
system embedded in a constant MF also possesses a constant of motion - 
"pseudomomentum". As a result for the case of zero total electric charge 
Q = the cm. motion can be removed from the total Hamiltoniai][^ The 
simplest example is a two-particle system with equal masses mi = 1112 = m 
and electric charges ei = —62 = e. We define 

R-=^^^Y^' ^ = ri-r2, P = Pi-f-p2. (1) 

^The case Q ^ is more complicated and will be considered elsewhere. 
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Straightforward calculation in the London gauge A = ^(B x r) yields 



One can verify that the following "pseudomomentum" operator F commutes 
with the Hamiltonian ^ 

F = P + |(Bx77). (3) 

This immediately leads to the following factorization of the wave function 
(WF) 

^(R, T7) = V7(r/) exp jzPR - i^(B x r7)R| . (4) 

For the oscillator-type potential V{ri) the problem reduces to a set of 
three oscillators, two of them are in a plane perpendicular to the magnetic 
field and their frequencies are degenerate, while the third one is connected 
solely with V{ri). 

Next we briefly elucidate the spin interaction in presence of MF. The 
corresponding part of the Hamiltonian may be written as 

Hs = 4:ahf {(71(72) - fiB{(Ti - (T2), (5) 

where ei = —62 = e > and fi > 0. Diagonalization of Hs yields the 
following four eigenvalues e.g. for uu system, comprising both p and vr levels. 



/ 



\ 



Tl , (6) 



where we assume that B is aligned along the positive z-axis and i? = |B|. In 
a strong MF when pB > cihf spin-spin interaction becomes unimportant and 
£^3^4 ^ ±2/ii?. For the lowest level E^^"^ this corresponds to a configuration 
I H — ) when the spin of negatively charged particle is aligned antiparallel to 
B, and the spin of the positively charged one - parallel to B. This means that 
the spin (and isospin) are no more good quantum numbers and eigenvalues 
^ correspond to the mixture of spin 1 and spin states. As a result the 
qq state will split into 4 states (two of them coinciding e[^^ = e!^^). Till 
now we treated a nonrelativistic system, to incorporate relativistic effects we 
shall exploit the path integral form of relativistic Green's functions [8]. 
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3 Relativistic qq Green's function and effec- 
tive Hamiltonian 



The starting point is the Fock-Feynmann-Schwinger (world-hne) represen- 
tation of the quark Green's function [8]. The role of the "time" parameter 
along the path zj];\si) of the i-th quark is played by the Fock-Schwinger 
proper time Si,i = 1,2. Consider a quark with a charge in a gluonic field 
and the electromagnetic vector potential A^^\ corresponding to a constant 
magnetic field Bi. Then the quark propagator in the Euclidean space-time 
is 

Siix, y) = {m, + d- igA - le.A^'^)-] = (m, + (7) 
The path-integral representation for Si is |8] 

POO ^ 

5,(x,y) = K-Z}«) / rfsi(Dz),,e-^»$«(x,i/) = K-D»)G,(x,y), (8) 



where 

<!>'^\x,y) = PAPpexp (^zg A.dzf + le, Af dzf^ x 

X exp (^j^ dTia^^{gF^y + eiB^u)^ . (10) 

Here F^y and B^^ are correspondingly gluon and MF tensors, Pa, Pf are 

) hold for the quark, 
i = 1, while for the antiquark one should reverse the signs of and g. In 
explicit form one writes 

<y,uF,. = ( aH ) ' "^^^^^^ = ( ^ ) • ^^^^ 

Next we consider qiq2 system born at the point x with the current Jri(2^) = 
qi{x)Viq2{x) and annihilated at the point y with the current Jr2(l/)- Here x 
and y denote the sets of initial and final coordinates of quark and antiquark. 
Using the nonabelian Stokes theorem and cluster expansion for the gluon 
field(see [TT] for reviews) and leaving the MF term intact, we can write 

G,,g,ix,y) = ds^j^ ds2{Dz^''^),y{Dz^^^),y{TWM))A^ 



ordering operators, a = 
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xexp(^ei r Afdz^;^-ie2 r A^'^ dzj^^ + f dn{aB)-e2 H dr^iaB)), 

Jy Jy Jo Jo 

(12) 

where 

f = tr{T^{m,- D,)T2im2- D2)), (13) 
and Fi = 7^, r2 = 'ju for vector currents, while 
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{WM))a = exp i^-^- J rf7r^,(l)rf7rA.(2)(F^.(l)FA.(2))j , (14) 

where dn^^ = ds^y + aj^Jdri — aj^JdT2, and dSfj_y is an area element of the 
minimal surface, which can be constructed using straight lines, connecting 
the points z''^\t) and z^\t) on the paths of qi and q2 at the same time 
t [U Unj. Then the spin-independent part of the exponent reduces to the 
confinement term Vconi{r) plus color Coulomb potential Vcoui, while spin- 
dependent part VsD depends also on proper time variables ti,T2, (see [H] 
for derivation and discussion). For the case of zero quark orbital momenta 
with the minimal surface, discussed above, one obtains a simple answer for 



{W„iA))A = exp i-j^JdtE 



_ -.(2)1 _ Z ^ 



0" ' — 



3|z(i) -z(2)| 



(15) 



containing Konf (|i"|) and Vcoui(|i"|)- Here a is the string tension, a = 0.2 GeV^ 
in our calculations. 

At this point we introduce the method of einbein variables (effective 
masses) defined via the connection between the proper time and the 
real time tf = Zi{Ti) 

dr, = / ds,{D^'h^%y = const J Du,{t){D^'^ z^%y. (16) 

In this way the path integral in Dz^^ is replaced by Dijj'^'^\ and the latter can 
be denoted as: / Dw^^^Dw^^^f...] = see [15] for the details. 

First we need to find the Hamiltonian Hq^q^ of the system at tf = tf = t^. 
To this end we define the Euclidean Lagrangian go- write — 



2uji-ij^ = 2uJiZk, k = 1,2,3. Then all terms in the exponents in (12), (14) 
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and (16) can be represented as exp(— / dt L^^^^) and thus we arrive at the 



following action 



S. 



Te 



+ E 



(ToB 



2uJi 2uj' 



-ei— \-e2- hcr|z^^^ 



.(2)1 



2ui 



3|z(i) 



Z(2)|2 

(17) 

Here A^^^ is the fc-th component of the QED vector potential, a is the 
QCD string tension. The next step is the transition to the Minkowski metric. 
This is easy since confinement is already expressed in terms of string tension. 
We have exp( 



H. 



J L dtE) — exp(z / L dtu), Hm, and 
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Pk 



dz 



(i) 



(18) 



Next comes the key point of the Method of Effective Masses P, HO]- It 
comprises the replacement of the path integral averaging over ui , U2 by the 
stationary point analysis. The applicability of this approximation may be 



justified by the following arguments. The qq Green's function (12) integrated 
over d^{x — y) takes the "heat-kernel" form 



Integrating (12) over d^{x — y), one obtains a simple expression: 



Gi',^lix,y)d'ix-y) 



2M, 



(A) 



(19) 



(20) 



Here e^^^ is the polarization vector for the polarization state A, and M'^\ 
are correspondingly the Hamiltonian eigenvalue and quark decay constant, 
2M^^ in the denominator stems from the normalization of the relativistic 
wave functions, n runs through all ordering numbers of the spectrum. All 
these quantities are functions of wi, uj2- Therefore the integral (20) may be 
symbohcally written as {Ke'^^'^)^^^^^ = J DuiDu2K{ui,U2)e-'^^'^'^''''^^^'^ and 
it is essentially defined by the region of the stationary point of the exponent. 

The effective masses cuj are to be found from the minimum of the total 
mass M{ui) , as it was suggested in [10]. To introduce the minimization 
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procedure and to check its accuracy we shall begin by the calculation of the 
eigenvalues of one and two quarks in MF, and the energy of the ground state 
of a relativistic charge in the atom in the next section, reproducing the known 
exact results. 

We have the following equations defining Ui from the total mass M{ui) 



dM{ui 
duji 



0. 



(21) 



For a single quark in MF the first of the above equations gives 

pl + ml + \eB\{2n + l)-eBa, u 
= + 2 • 

Then the second equation yields the correct answer 

M„ = {pI + ml + \eB\{2n + 1) - eBa^f'^. 



(22) 



(23) 



Now we turn to the case of qiq2 system and introduce the coordinates which 
are the generalization of ([T]) 



R 



UJ1 + UJ2 



z«-z(2), 
dr]' 



(24) 



(25) 



It is convenient to introduce the following two additional parameters 



U1U2 UJi — OJ2 

s 



UJI + UJ2 



UJi + UJ2 



(26) 



As before, for simplicity we consider only the neutral meson, so that ei 
—62 = e. Then the total Hamiltonian may be written as 



H, 
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(27) 



where 



2uJi 



— P + TT- -B X (R 

UJ2 2 



UJ 



n 2 



+ 



77 + -B X (R 



UJ 



UJ2 



77 



1 

2^ 



-B X R + s-B X 77 
2 2 ' 



+ 



1 



2(ui + U2) 



P - -B X ry 
2 ' 



(28) 



Equation (28) is an obvious generalization of ([2]). The two other terms in 
1^ read 

-"cr = \ , y^y) 



W 



2ui 

Konf + V^Coul + AVT 



a\r]\ 



2U2 
3 T] 



+ AW, 



(30) 



and AW contains self-energy and spin-spin contributions. One can verify 
that the "pseudomomentum" operator in ^ introduced in Section 2 com- 
mutes with Hb and hence we can again separate the cm. motion according 
to the ansatz Q. Then the problem reduces to the eigenvalue problem for 
ip{r]) with the Hamiltonian Hb having the following form: 



1 f d e„ ^ 



2(ui + UJ2] 



(P - eB X -qf 



(31) 



For P X B = the system has a rotational symmetry and the cm. is freely 
moving along the ^-axis. Here we shall consider a state with zero orbital 



momentum (L^)^ = 77 x 
internal space operator 



idT] 
1 

" 2Q 



0. As a result Hb is replaced by a purely 
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To test our method we put W = and arrive at the equation 

(Ho + K)ip = M{uJi,U2)^p. 



(32) 



(33) 



Consequent minimization of M{ui,U2) in uji,U2 , as in (23) , yields the 
expected answer for the two independent quarks, 



M 



ml + eB{2ni + 1) - ecriB + Jm^ + eB{2n2 + 1) + ecraB. (34) 



Treating confinement and color Coulomb 
terms. The magnetic QCD collapse 



From(30), (32) it is clear, that inclusion of Konf and Vcoui iia Hq + W leads to 



a differential equation in variables f]±,r]z, which can be solved numerically. 
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However, in order to obtain a clear physical picture, we shall represent Konf 
in a quadratic form. This will allow to get an exact analytic solution in 
terms of oscillator functions with eigenvalue accuracy of the order of 5%. 
The color Coulomb contribution will be estimated as an average Vcouilv^), 
thus yielding an upper limit for the total mass. 
For V^onf we choose the form 



X^conf ^ Konf = o - + 7 (35) 




Here 7 is a positive variational parameter; minimizing Konf w.r.t. 7, one 
returns to Konf- We shall determine M{uiU2,'y) corresponding to Konf, and 
to define 7 an additional condition 



dM{uji,uj2,'j) 



(36) 



7=70 



will be added to (21). As a result M(co'i'^\ Ci;2°\ 7o) will be the final answer 
for the mass of the system. The difference of the exact numerical solution 
from that obtained with the genuine potential Konf does not exceed 5% . 
The solution of the equation {Hq + V^onf)v' = M{ui,U2,'yo)v fo^ the ground 
state is 

where r± = ^J~^ (l + ^^2^2 ) ^ , 1^0 = (^^^ ^ . As we shall see below, for 
the lowest mass eigenvalue with eB ^ a, one has r± ^ ^7=, ro ~ ^ and 
the (gi, ^2) system acquires the form of an elongated ellipsoid. Similar quasi- 
one-dimensional picture was observed before for the hydrogen-like atoms in 
strong MF (SJ E] ■ In such geometrical configuration Vcoui manifests itself in 
a peculiar way, again similar to what happens in hydrogen, or positronium 
atoms. However, as we shall see now, in QCD, at least in the large Nc limit, 
the outcome is drastic. 

We turn now to the color Coulomb term. As a starting point we present 
another check of our approach, namely we shall obtain the ground state 
energy of two relativistic particles with opposite charges without MF inter- 
acting via the Coulomb potential. The corresponding Hamiltonian reads 
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H = Ho + - then H(f) = M(j), and for eS = we have 

M = -^ + ^^i±^ + ^^i±M. (38) 
2 2wi 2^2 

Minimizing in Ui in the hmit m2 ^ mi (the hydrogen atom), one obtains 

M = miVl - a2, (39) 

which coincides with the known eigenvalue of the Dirac equation. 

In our (91^2) case one can calculate the expectation value of Vcoui = 
— I with the asymptotic freedom and IR saturation behaviour in p-space 



see 



TQI for a short review) 



4:71 

aJq) = (40) 



(f{rl - rl) 



where Mb is proportional to y/a, Mb ^ 1 GeV (TB] . With the wavefunction 
(37) the average value of Vcoui takes form 

Vconi{q)^\q)j;^ = -7^ / as{q)dqe- — I 
[Zir)-^ Sir Jo 

where /(a^) = J^^ dxe " ^ . Estimating the integral in (41), for eB ^ a, i.e. 
for To ^ r_|_ one obtains for massless quarks 

AMcoui ~ -i^iriln^ ^ -v^lnln — . (42) 

6r oPo 'fj_ a 

With eB increasing the upper bound for the qq mass is boundlessly de- 
creasing. The exact eigenvalue should lie even lower. Surmising (as will be 
confirmed in the next section) that the contribution of the remaining part 
of Hamiltonian to the total mass is Mq 2-(/o", we can estimate the upper 
limit of eB compatible with the conditions Mq + AMcoui ^ namely 

ieB)2!^^ ^ aexp (exp ( J J ~ 2.5 ■ lO'^^G ^ 2.8 ■ lOV. (43) 



As we shall see, numerical calculations yield much smaller limit: {eB)'^^ ^ 
6 GeV^. We note, that this upper limit is much smaller, than obtained in 
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QED from the positronium collapse [S]: 



(e5)Sfx^ = ^ exp + 2Ce) - 10^°G. (44) 



One should note, that the QCD limit (43) is unaffected by higher order 
gluon loop corrections, since those contain only gluons, not sensitive to MF. 
However, the quark loop corrections to Vcoui are growing like \eB\ and can 
possibly ensure the necessary screening. This is in line with QED, where such 
corrections produce screening and stabilization in the hydrogen case [6l [7] . 



Therefore the magnetic QCD collapse (42) refers only to the large Nc limit, 
when the quark loops contribution can be disregarded. 

We shall not elaborate here more on this problem and its significance, 
leaving the topic to a dedicated paper. 



5 Meson masses in magnetic field 

Our next task is to calculate analytically the mass M„(u;i, u;2, 7) of a (^1^2) 
meson. We have to solve the equation 



(Ho + H^ + W)^n{v) = M„(a;i,a;2,7)^n(r/), 



(45) 



where Hq, H^, W are given in (29 32). The result is 

Mn{Ui,U2,j) = £nx,n. + ^ ^ ^ ^- + {A Mconl) +AMsE , 



2Wi 



2cjo 



(46) 



where 



1 

2^ 



W + :^(2., + i) + 

7 V 7 



1 

riz -\ — 
2 



+ 



70" 



(47) 



AMcoui is given by (41), while according to [HI [T2] Vss and AMse are 
given by 



V, 



ss 



8tt ahfl^Pn 



9 U1 + U2 



{(Ti(T2),AMsE 



2a 



2a 



(0) (0) • 

TTLOl 7TUI2 



(48) 



We note that both V^^ and AMse are to be considered as corrections 
and contain Ui^\uj2'\ obtained from minimization of the remaining part of 
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the Hamiltonian. The parameter 70 (see (35)) is defined from the condition 

dM 



7=70 



9e 
d'y 



0, 



(49) 



7=70 



In the lowest state \u l,u t) we have ctiB = — cr2B = B, 



(0) 



0;^°^ and w^"^ is obtained from 



dM 



^2 



(0) 



0, i = 1,2. 




Figure 1: The mass of the system as a function of y/eB. See the text for 
explanations. 

In Fig.l we plot the mass of the system as a function of y/eB. Calculations 
were performed according to (46) and the minimization procedure. The solid 
curve corresponds to the configuration \u l,u^), the dashed one — \u^,u^). 
The circles are from lattice calculations [T7j, the squares — from [TS]. One 
can see that the mass is slowly decreasing and reaches zero at (e-B)^^f (note, 
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that the results plotted in Fig.l were obtained for massless quarks). The 
behaviour is in agreement with the "magnetic QCD collapse" phenomena 
discussed above. 



6 Discussion and conclusions 

In our treatment of relativistic quark-aniquark system embedded in MF we 
relied on pseudomomentum factorization of the wave function and relativistic 
einbein technique. The Hamiltonian for neutral mesons in MF, containing 
confinement, colour Coulomb and spin interaction was derived. Using a suit- 
able approximation for confining force we were able to calculate analytically 
meson masses as functions of the MF. Our eye was predominately on the low- 
est level with its mass decreasing with MF growing. This state is a mixture 
of the p° and 7r° as can be seen from its spin and isospin structure. Indeed, 
uu system under consideration is a mixture of isospin 1 = and 7 = 1 states, 
and at large MF it has a spin structure \u '\,u 1), which is a mixture of 5* = 
and S = 1 states. With MF growing the mass of this state decreases (see Fig. 
1), while the masses of all other states increase as \/eB. A significant point 
is that these results are in line with recent lattice simulations [T71 [18] (see 
Fig.l). Calculating the Coulomb energy (AMcoui) we have obtained the un- 
bounded negative contribution at large MF proportional to ( — y^lnln — 



which makes the total mass neg ative for eB > e^^S^ ^ 10 GeV^ 

Unlike the situation in hydrogen-like atoms, where loop corrections are 
able to produce saturation E], in QCD gluon loops are MF blind and 
uncapable to improve the results, while quark loops are suppressed in large Nc 
limit. We call this problem "magnetic collapse in QCD" and plan to discuss 
it in detail in a separate paper where quark loops and gluon polarization 
operator will be considered, and possibly improve the situation, similarly to 
the hydrogene atom case. 

In this paper to simplify things we started with p° meson states at 



-8 = taking 7^ in place of Fi and r2 in (13). In this way we essentially 
left aside the complicated problem of chiral dynamics and pseudo-Goldstone 
spectrum. As explained above in this oversimplified picture we can consider 
the lowest state as a mixture of p° and 7r°. These two states are splitted by 
hyperfine interaction, this splitting is insignificant in the large MF limit. It 



^In fact uu is a mock p° 
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is legitimate to compare the results of such treatment with the lattice data 
pn [18] since in the latter the quark masses are not small and thus the chiral 
facet of 7r° is suppressed. As shown in Fig.l our analytical results are in 
agreement with lattice calculations [T71 [18] both for p"{uu) states | ti> and 
I tt>- The behavior of the total mass Mq supports the conjecture of the 
"magnetic QCD collapse" existing in absence of quark loop corrections. 

The methods used above can be generalized to the charged states thus 
shading the new light on the problem of charged vector boson condensation 
suggested in [7]. As a preliminary foresight we note, that instability is an 
inherent property of elementary spin 1 bosons, while p-meson can not be 
considered as such an object, when MF is so strong that the Larmour radius 
is equal or smaller than its size. Another system which can be treated using 
the same technique is the neutral 3-body system, like neutron. The results 
might be important for the neutron stars physics. 

The authors are grateful for useful discussions with V.A.Novikov, M.I.Vy- 
sotsky and S.I.Godunov. We are indebted to V.S.Popov for important re- 
marks. B.K. gratefully acknowledge the support RFBR grant 10-02093111- 
NTSNIL-a. We are pleased to thank M.Chernodub for his remarks and 
suggestions in response to v.l of this paper. 



References 

[1] Proceedings of the international Conference on Physics in Intense 
Fields (PIF 2010), November 24-26, 2010 KEK, Tsukuba, Japan, 



http://atfweb kek.jp/pif 2010/. 



[2] J.M.Lattimer and M.Prakash, Phys. Kept. 442, 109 (2007). 

[3] D.E.Kharzeev, L.D.McLerran and H. J.Warringa, Nucl. Phys. A803, 227 

(2008) ; 

V.Skokov, A.Illarionov and V.Toneev, Int. J. Mod. Phys. A24, 5925 

(2009) . 

[4] T.Tajima, Eur. Phys. J. D55, 519 (2009). 

[5] A.E.Shabad and V.V.Usov, Phys. Rev. Lett. 98, 180403 (2007), Phys. 
Rev. D73; 125021 (2006). 



15 



[6] B.Machet and M.I.Vysotsky, Phys. Rev. D83: 025022 (2011); 
S.I.Godunov, B.Machet and M.I.Vysotsky, |arXiv:1112T8 91 [hep-ph]. 

[7] J.Ambjorn and P.Olesen, Nucl. Phys. B315, 606 (1989); Phys. Lett. B 
218, 67 (1989); 

M.N.Chernodub, Phys. Rev. Lett. 106: 142003 (2011), Phys. Rev. D82: 
085011 (2010). 

[8] Yu.A.Simonov, Nucl Phys. B307, 512 (1988); Yu.A.Simonov and 
J.A.Tjon, Ann. Phys. (N.Y.) 300, 54 (2002). 

[9] L.Brink, P.Di Vecchia and P.Howe, Nucl Phys. B 118, 76 (1977); 
Yu.S.Kalashnikova and A.V.Nefediev, Phys. At. Nucl. 61, 785 (1998). 

[10] Yu.A.Simonov, Phys. Lett. B 226, 151 (1989); A.Yu.Dubin, 
A.B.Kaidalov and Yu.A.Simonov, Phys. Lett. B323,41 (1994); Phys. 



At. Nucl 56, 1745 (1993); hep-ph/9311344 



[11] A.Giacomo, H.G.Dosch, V.LShevchenko and Yu.A.Simonov, Phys 
Re pt. 372, 319 (2002); Yu.A.Simonov, Phys. Usp. 39, 313 (1996) 
|arXiv:hep-ph79709344 ; Yu.A. Simonov, QCD and Theory of Hadrons in 
"QCD: Perturbative or Nonperturbative." Interscience, Singapore, 2000 
rhep-ph/9911237, 

[12] A.M.Badahan and B.L.G.Bakker, Phys. Rev. D66, 034025 (2002); 
A.M.Badahan, B.L.G.Bakker, Yu.A.Simonov Phys. Rev. D66, 034026 
(2002); A.M.Badahan, B.L.G.Bakker, Phys. Rev. D84, 034006 (2011); 
A.M.Badahan and B.L.G.Bakker, Phys. Rev. D81, 071502 (2010). 

[13] W. E. Lamb, Phys. Rev. 85, 259 (1952); L.P.Gor'kov and LE. 
Dzyaloshinskii, Soviet Physics JETP, 26, 449 (1968); 
J.E.Avron, LW.Herbst, and B.Simon, Ann.Phys. (NY), 114, 431 (1978); 
H.Grotsch and R.A.Hegstrom, Phys. Rev. A4, 59 (1971); 

[14] Yu.A.Simonov, NucL Phys. B324, 67 (1989); A.M.Badahan, 
A.V.Nefediev and Yu.A.Simonov, Phys. Rev. D78, 114020 (2008). 

[15] A.M.Badahan, B.L.G.Bakker, Yu.A.Simonov Phys. Rev. D 75, 116001 
(2007). 



16 



[16] Yu.A.Simonov, Phys. At. Nucl. 74, 1223 (2011); [i^^v: 1011 .53861 [hep- 
ph]. 



[17] E.V.Luschevskaya and O.V.Larina, arXiv: 1203.5699 [hep-lat]. 
[18] Y.Hidaka and A.Yamamoto, arXiv: 1209.0007 [hep-lat]. 



17 



